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Abstract 

In practice daily volatility of portfolio returns is transformed to longer 
holding periods by multiplying by the square-root of time which assumes 
that returns are not serially correlated. Under this assumption this pro- 
cedure of scaling can also be applied to contributions to volatility of the 
assets in the portfolio. Close prices are often used to calculate the profit 
and loss of a portfolio. Trading at exchanges located in distant time zones 
this can lead to significant serial cross-correlations of the closing-time re- 
turns of the assets in the portfolio. These serial correlations cause the 
square-root-of-time rule to fail. Moreover volatility contributions in this 
setting turn out to be misleading due to non-synchronous correlations. We 
address this issue and provide alternative procedures for scaling volatility 
and calculating risk contributions for arbitrary holding periods. 
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1 Introduction and Motivation 

In this article we consider a portfolio consisting of n > assets whose returns at 
time t £ Z are modelled by the random vector X t := (X*, . . . , X^) T for t G Z, 
where Z denotes the integers. We will use methods from time series analysis 
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where Z is the natural range for the time index (see for example Mills 1993 ). 
The percentage weights of the assets in the portfolio are denoted by the vector 



Thus we allow for short-selling as well as leverage 



where the leverage can be measured by J2i=i ^i- 

Building a portfolio using the weights A the random portfolio return at time 
t € Z, denoted by X t (X), is given by 



X t (X) := \ T X t = \Xl for t € Z. 



(1.1) 



For our analysis we will assume that {Xt)te.Z and therefore, as we show in 
Proposition 2.6 also (X t (X))t<£Z is a weakly stationary multivariate respectively 
univariate time series. Thus mean and covariances exist and covariances depend 



on the lag and not the absolute time points (for details see for example Box et al. 

' 1986] ) 



2008 , Brockwcll and Davis 



In practice it is assumed that the component time series in {X t )t£i do not 
exhibit any serial correlation. Thus it is usually assumed that yesterday's return 
of asset i is not correlated with today's return of asset j. But there are situations 
where this assumption is clearly not fulfilled as we will sec in the following. 

Holding a globally diversified portfolio of stock-index futures and analyzing 
the returns between close prices one sometimes experiences significant correla- 
tions in lagged returns. This may be due to non-overlapping trading times of 
various exchanges (e.g. New York and Tokyo). For a detailed analysis of the 
correlation bias we refer to Kahya 1997 , Coleman 2007 . This problem can 



not be addressed by taking prices from a point in time where all exchanges in- 
volved are open as such a moment might simply not exist. This lead-lag effect 
is a prime example for the application of multivariate time series analysis, see 



for example |McNeil et al. 2005 Chapter 4, Section 5] and deJong and Nijman 



1997 for empirical aspects in high frequency. 



Modelling prices in different time zones we consider days t G Z and closing- 
time fractions x\ < xi < • • • < x n of the n exchanges where the assets in our 
portfolio are traded. By this we mean that the first market closes each day at 
the point in time x\ (e.g. if x\ = 1/3 then we mean 08:00 CET) and so forth. 
For the return series of a portfolio (X t (A)) tG z we have to define the following 
notiorQ 

Definition 1.1 (Closing-time return). Calculating the return of an investment 
using ( 1.1 ) where each return X\, i = 1, . . . , n is calculated from the close price 



in the respective market PI, i.e. 



X\ = {P i 



t+Xi 



PI 



1 + Xi, 



for t e z, 



(1.2) 



where Xi is the closing-time fraction, is called closing-time return. 



Thus we use the notion of closing-time return as the return on investment 
that comes from book-keeping the prices each day at the point in time when 
the respective market closes. This is e.g. 08:00 a.m. CET for Japanese stocks 
and 10:00 p.m. CET for stocks traded in the US. This return is often used by 



1 We thank an anonymous referee for pointing out that the following clarification of notions 
is essential. 
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accountants to calculate the performance of an investment fund or the profit 
and loss of a trader. 

As one can see from the example of Japan and USA, a closing-time return 
does not include the reaction of markets to events that happen when they are 
closed. Economically an event that affects the US market while the Japanese 
market is closed will potentially affect Japan one day later. Profits and losses 



accounted for by (1.2) affect the aggregate performance of the portfolio with 
some lag as profits and losses can be increased or decreased by the reaction of 
the Japanese market on the following day. 

This lagging reaction changes our view of the risk that we economically have. 
The idealized concept of a return that takes into account lagging reactions on 
one and the same day is called contemporaneous return in this article as opposed 
to the non-contemporaneous closing-time return. This return could be thought 
of as a result of shifting all trading times, i.e. opening and closing-times, such 
that trading takes place contemporaneously. Note that the trading times in 
the US and Japan do not even overlap. Thus the calculation of the risk in 
contemporaneous returns can be seen as an estimation problem that can for 
example be tackled by the Newey-West estimator (see Newey and West |1987 ) . 
We briefly show an application of this estimator in the end of Example [1] below. 

Summing up we have a closing-time return used by accountants to asses the 
performance of a fund or a trader in a straight and clear way. On the other 
hand we have the contemporaneous return that tries to reflect the economic in- 
teraction of markets. In this paper we mainly deal with closing-time returns and 
provide formulas to calculate genuine risk and risk contributions. In the follow- 
ing we give a concrete example and motivate our focus on vector-autoregressive 
models and especially the VMA(l) case. 

Example 1 . Consider a portfolio of stock index futures traded in Japan, China, 
Europe, South Africa and the United States. Note that weights sum up to 105% 
which corresponds to a small degree of leverage. The following table shows the 
names of stock indices in these regions, example weights and the closing-times 
of the exchanges in CET. For such a portfolio there is no point in time where 
all exchanges are open. Using close prices at all exchanges one has to treat 
the lagged correlations between the returns of different assets as it is illustrated 
below. 



Asset 


Currency 


Exposure 


Closing (CET) 


1 


Topix (TOP) 


JPY 


15% 


08:00 a.m. 


2 


H-shares (HSH) 


HKD 


15% 


10:00 a.m. 


3 


DJ Euro Stoxx 50 (DJS) 


EUR 


15% 


05:30 p.m. 


4 


Swiss Market (SWI) 


CHF 


15% 


05:30 p.m. 


5 


JSE TOP 40 (JSE) 


ZAR 


15% 


06:30 p.m. 


6 


Russell 2000 (RUS) 


USD 


15% 


10:00 p.m. 


7 


NASDAQ 100 (NAS) 


USD 


15% 


10:00 p.m. 



The correlations of the returnsjj of the indices underlying these futures con- 
tracts are plotted in Figure[l]where the numbering corresponds to the numbering 
in the table above. At lag zero we see that geographically close markets have 

2 Estimatcd from a 250 days window ending in April 2010. 
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relatively strong correlation in returns. One recognizes an Asian (TOP,HSH), a 
European (DJS,SWI) (also showing significant correlations with South Africa) 
and an American block (RUS,NAS). In Figure [I] the correlations of the returns 
(y-axis, today) to lagged returns (x-axis, yesterday) are plotted. The auto- 
correlations on the main diagonal (i.e. correlations of returns of one asset with 
lagged returns of the same asset) are not significant which is consistent with a 
random walk. But on the other hand this plot illustrates that, in our exam- 
ple, returns in Asian markets have significant correlation with lagged returns of 
the European as well as the US block. This illustrates the lead-lag effect. We 
will come back to this portfolio in Section [3] and apply advanced procedures to 
analyze risk. 

We continue with a rather heuristic analysis of the situation in the spirit of 



Kahya 1997 , Coleman 2007 and Bergomi 2010 . Let 



(B t ) tf zR - (B t , B") t( z R 

be an n— dimensional Brownian motion with the identity matrix as correlation 
matrix, i.e. 

Cov(Bi,Bi) = 0, if Mi. 

Then, with S = (ot,j)?j=i an n x n positive-semidefinite matrix and E 1 / 2 its 
Cholesky-decomposition, it is clear that (Btjtem. defined by B t = E 1 / 2 • B t for 
t € M is a Brownian motion with covariance matrix E, meaning that 

Cov(B l t ,B{) =ta itj , for t > and i,j = l,...,n. (1.3) 

Then the increments B t — B s have a multivariate normal distribution with co- 
variance matrix (t — s)E. Recording prices at different closing-times during 
the day can be described by observing prices at integer times t shifted by the 
closing-time fraction Thus we model our closing-time returns 

from Definition [Ll] by the vector 

x t = {X t , A") = (B t+Xl - Bl +Xl _ 1: B" +Xn - B? +Xn _ 1 ), for t G Z, 

where we assume that assets are ordered with respect to their closing-times, 
i.e. x\ < xi < ■ ■ ■ < x n . If x\ = Xi ■ ■ ■ — x n then we can shift time by x\ by 
stationarity of Brownian motion and get back the standard model with E as 
the covariance matrix of the asset returns. Assuming without loss of generality 
that Xi > Xj, i.e. that market i closes later than market j then we get for the 
general case 

Cov(Xi,Xi)=Cov(Bi +Xi -Bi_ 1+Xi ,Bi +x . -Bl 1+X .) 
= Cij(l - (%i - Xj)), 

where we have used 

\t — 1 + Xj , t + Xi] = [t — 1 + Xj , t — 1 + Xi) U [t — 1 + Xi, t + Xj) U [t + Xj,t + Xi] 
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as a decomposition for the time interval and 

B l+ Xl - Bl-x+xi = B l+ Xj - B l-i+ Xl + B l+ Xz - B l+ xp respectively, 

S v " 

B\ on [t — l-\- X i,t-\-Xj) 
— R- 7 — R J — R J -l-fi- 5 — R- 7 

" V ' 

B{ on [t-l+Xi,t+Xj) 

for the Brownian motion. Thus the covariance matrix of the closing-time returns 
X t is given by 

:=Ki(l-|*i-Zil))" j=1 - (1-4) 

Note that on the main diagonal we have Oi.i = of, i.e. variances are not biased 
but covariances are. 

Furthermore we consider the lag one covariance matrix of X t . Let 

r(i) = Cov(x t+1 ,x t ) 

such that 

r(l) 4J = Cov(Xi +1 ,X(), for i,j ! 1 1 . 

Then, assuming Xj > Xi , due to the overlap of the interval [t — 1 + Xj,t + Xj] 
where we observe X\ , respectively [t + Xi,t + 1 + a^] where we observe X| , lt we 
have 

m cv< yJ\ J ' if Xj < x i? and, 

Cou(X 4 % 1; X/) = i (1.5) 

With the assumed ordering of the markets T(l) is an upper triangular matrix 
with zero main diagonal. The covariance matrices of lag two and greater are 
zero as there is no overlap in the time intervals of interest anymore. For example 
at lag two consider that 

Cov{Xl +2 ,Xi) = Cov(Bl +2+Xi - Bi +1+Xi ,B{ +Xj - B\ +X] _ x ) = 0, 



as 



[t - 1 + Xj,t + xj] n [t + 1 + Xi,t + 2 + Xi] = 



Finally note that the vector (X t +i,X t ) is jointly Gaussian and the conditional 
law of Xt + \ given X t is Gaussian with mean T(l)Yi~j}X t and covariance matrix 

s x - r(i)s^r(i) T . 

Thus knowing the conditional mean of X t+ i given X t we can write 

X t+1 = $X t + £ t+1 , forteZ, (1.6) 

with $ = r(l)!]- 1 and a process {£t)t& which is not a white noise process as 
some lines of calculations show. Although the (visual) form (1.6 1 reminds of a 
VAR(l) model, the vanishing lag covariances greater than one are inconsistent 
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with a VAR(l)-model (see the form of the lagged covariance matrices in (3.4) 
in Section [3]). We rather assume a VMA(l)-model of the form 



%t+i — ©-£■( + Zt-\ 



for t G Z, 



(1.7) 



with a white noise process (Z t )t^L- We will come back to this model class in 
Section [3l 

As mentioned before calculating the covariance matrix of contemporaneous 
returns leads to a special estimation technique and we give two examples. It is 
straight forward to define an estimator by 



+r(i) + r(i) J 

:=r(0) 



(1.8) 



It is known that the estimator ( 1.8 ) does not always yield a positive-semidefmite 
covariance matrix (see Newey and West 1987| ). However, considering Equa- 
tions (1.4 1 and (1.5) we see that this covariance estimator of contemporaneous 



returns, E, equals X asymptotically, which is the covariance matrix of the Brow- 



nian motion defined in (1.3). Thus in this setting the estimator (1.8) reveals 



the true contemporaneous covariance matrix. For the general case there is an 
estimator which always yields a positve-semidefinite covariance matrix, the so 
called Newey- West estimator. The estimator introduced in Newey and West 
1987 up to lag 1, denoted by £jvw> is given by 



W-r(o) + ^(r(i) + r(i) T ) 



(1.9) 



However in the setting of non-contemporaneous trading we prefer (1.8) 




Figure 1: Heatmap of correlations of lag zero. 

In the following we mainly consider closing-time returns and do not always 
mention this fact. If we make some remarks concerning contemporaneous re- 
turns this will be mentioned explicitly. 
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Figure 2: Heatmap of correlations of lag one. 



The risk of the return X t (A) holding the portfolio for one day is often mea- 
sured by volatility which is simply the square-root of variance. Stressing that 
this volatility depends on the asset weights A, we write 



a(X) := y/Var(X t (\)) 



\ 



Var(£>Xj) 



for t € Z. 



(1.10) 



Remark 1.2. If we furthermore assume that the vector X t is normally distributed 
for each t, then knowing the volatility and the mean we can calculate value-at- 



risk and expected shortfall by the well-known formulas McNeil et al. 
Chapter 2]. 

Remark 1.3. If £ 



2005 



denotes the covariance matrix of X t 



t, i-e. 



Cov(Xl, X 3 t ), then it is well known that 



<r(A) = VFea. 

Regulatory rules require the calculation of the risk of a portfolio for a holding 
period of d > 1 days (e.g. 10 or 20 days) and it is common to quote volatility 
as per annum (i.e. d — 250 or d — 252). For the clearness of presentation we 
will assume that we model returns at t — 1, . . . ,d when considering a holding 
period of d days. 

Definition 1.4 (Volatility for a holding period of d days). We denote the 
volatility of the portfolio return given in (1.1 1 holding the assets for d days by 
c(A, d) and thus 



(x(X,d) 



\ 



War(Y^X k (X)). 



(1.11) 



k=l 



In practice it is assumed that the returns X t (X) are serially uncorrelated for 
t E Z and that variance is stationary, i.e. that War(Xk) — c(A) for k = 1, . . . , d. 
These assumptions lead to the square-root-of-time rule. 

Proposition 1.5 (Square-root-of-time rule). Under the above assumptions it 
holds that 

a(X,d) = a(X)Vd, ford>l, (1.12) 
where cr(A) denotes the one-day volatility. 

Note that, due to its simplicity, the square-root-of-time rule is often used 
in practice. Using it without checking its appropriateness can lead to poor risk 
estimation (see McNeil et al. , 2005| Chapter 2] and references therein) . 

Having calculated the volatility of the portfolio return the Euler allocation 
is used in practice to define risk contributions by the assets (see Tasche |2000 



2008 and McNeil et al. 2005[ Chapter 6] 



Definition 1.6 (Contribution to volatility). The contribution to the volatility 
of the portfolio return X t (X) of asset i by the Euler rule, denoted by Ci(A), is 
given by 

o-i(\):=\i for i = 1, . . . , 7i. (1.13) 

v /Var(X t (A)) 

Remark 1.7. Using the same notation as in Remark|1.3|one can easily see that 



Pi (A) = for i = 1, ...,n, 

cr(A) 

where (£A)j denotes component i of the vector SA. 

Assuming that asset returns X t are not serially correlated for t 6 Z we can 
make the following definition for the risk contribution for a holding period of d 
days, which we denote by <7j(A, d) for i = 1, . . . , n: 

Proposition 1.8 (Square-root-of-time rule for risk contributions). Under the 
above assumptions it holds that 

<7j(A, d) = <7j(A)Vrf, /or i = 1, • • • , n and d > 1, (1-14) 



where CTi(A) is given in (1.13). 

Remark 1.9. Defining cr(A) and o"j(A), i = 1, . . . , n as above it is easily seen that 
we have full allocation for all holding periods d, i.e. 

n 

Q-(A, d) = o~i(A, d), for d > 1. 
t=i 

Furthermore considering relative risk contributions it is easily seen by the follow- 
ing equation that they do not change when considering longer holding periods: 

Qi(A, d) _ <Ji(\)Vd _ gj(A) 
cr(A,d) " cr(A)\/d ~ o-(A) ' 

for t = 1 , . . . , n and (i > 1 . 
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The above rules for scaling volatility are well-known and used in practice. 
In the following section we drop the assumption that asset returns are serially 
uncorrelated. In Section [2] we derive general formulas for scaling volatility and 
volatility contributions in a setting with serial correlations. Modelling the port- 
folio return as a univariate process with auto-correlations we propose proper 
volatility scaling in Subsection |2.1| while the multivariate setting of Subsec- 
tion |2.2| additionally allows for the calculation of proper risk contributions. In 
Section [3] we come back to the concrete problem of non-contemporaneous trad- 
ing and apply these results to VMA(l)-models and derive the corresponding 
formulas. Then we return to the data of Example [T] and show that the square- 
root-of-time rule can seriously underestimate volatility. Furthermore it can give 
misleading indications of risk contributions. Finally we take a short detour and 
propose VAR(l)-models to tackle the problem of genuine auto-correlations in 



the sense of Anderson et al. 2005 



While we focus on auto-regressive modelling of the returns, a GARCH- 



approach to the problem of scaling volatility is given in Diebold et al. 1997 



and the scaling in a model with jumps in returns is considered in Danielsson| 
and Zigrand 2006 . For a study on the square-root-of-time rule for tail risk 
we refer to the recent paper Wang et al. 2010 . While an analysis of temporal 



aggregation of ARMA-models with GARCH errors is given in |Drost and Nijman| 



1993 we focus on risk contributions. To our knowledge the question of scaling 



volatility contributions is so far not covered in the literature. 



2 Volatility Contributions under Serial Correla- 
tion 

Recall that we analyze volatility scaling and volatility contributions in the set- 
ting of weakly stationary processes. There are circumstances where the assets 
that constitute the portfolio are not known and one can only model volatility on 
the level of portfolio returns. In this setting we will propose a scaling rule that 
takes auto-correlations into account. Later on we analyze the situation when 
assets are known. Then a bottom-up modelling of the portfolio is possible and 
the whole covariance structure of lagged asset returns can be estimated. 

2.1 The univariate model with auto-correlations 

First we make the following definition for the auto-covariance and auto-correlation 
function of portfolio returns (X t (X)) t £%: 

Definition 2.1 (Auto-covariances of a univariate time series). Let (X t (\)) tG % 
denote a univariate weakly stationary stochastic process, then the auto-covariance 
function and the auto-correlation function is denoted by 

j(k) :=Cov(X t (X),X s (X)) and (2.1) 
p(k) Cor(X t (X),X s (X)) = j(k)/ 7 (0), (2.2) 

for t, s £ Z where \t — s\ = k, respectively. 

In the presence of auto-correlations in the portfolio returns the scaling by 
the square-root of time is not accurate and the following proposition states the 
correct scaling in this setting. 
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Proposition 2.2 (Volatility for a holding period of d days). Let (Xt(X))tez be 
a univariate weakly stationary stochastic process with auto-covariance function 



j(-) as defined in Definition 2.1 then the volatility of the return when holding 
the portfolio over d > 2 days is given by 



a(\,d)=a(£x k (X)) 



fc=i 



\ 



ci-l 



d 1 (0) + 2j2(d-kh(k). 



(2.3) 



fc=i 



The proof of Proposition 2.2 is given in the appendix and the correction to 



scaling by the square-root of time is clearly seen in the following corollary. 

Corollary 2.3 (Scaling rule for the univariate model). Let (X t (X))tei, as in 
Proposition 2.2 and p(-) be its auto- correlation function then 



a (X,d) = cr(X)S(d), 



(2.4) 



where cr(A) is the one-day volatility from ( 1.10[) and the factor 6(d) is given by 



S(d) = 



\ 



d-1 



d + 2j2( d - k )p(k), ford>2. 



(2.5) 



k=l 



The short proof of this corollary can also be found in the appendix. Con- 



sidering (2.5) we see that in our univariate model the scaling is given by the 



square-root of time corrected by an expression taking into account all relevant 



auto-correlations. If these are zero then (2.5) reduces to the well-known for- 



mula (1.12) 



The next corollary gives a crude estimate of how much higher the true scaling 
factor can be compared to the square-root-of-time rule. 

Corollary 2.4 (Error when using the square-root-of-time rule). Let (X t (X))t^z 
as in Proposition 2.2 and S(d) be the correct scaling factor, then the proportion of 
correct scaling to an application of the square-root-of-time rule can be estimated 
as 



y/d 



< Vd, for d>2. 



(2.6) 



Proof. Considering (2.5) and noting that \p{k)\ < 1 for all k G Z, we get 



5(d) < 



\ 



d+2j2(d-k) =d, 



fc=i 



where we apply that ^2t=\(d ~ k) 



W 2 



d). 



□ 



The above corollary states that the proportion of the correct scaling factor 
and the square-root of time when the necessary assumptions are not fulfilled 
grows with the square-root of time. This conservative estimate does not assume 
anything non-trivial about the auto-correlations. In Figure [3] in Section [3] we 
will see that the concrete picture is not always that bad. 
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2.2 The multivariate model with auto-correlations 

Next we will analyze the situation if the constituent assets are known. In this 
setting a bottom-up modelling of the portfolio structure is possible. For the the 
covariance structure of the asset returns (X t ) te % we need the following definition: 

Definition 2.5 (Covariance matrix function of a multivariate time series). Let 
(Xtjtez denote a weakly stationary process in W 1 . Then we denote the matrices 
of serial covariances of lag fc = 0,l,2,...by 

r(fe) :=Cov(X t+k ,X t ), (2.7) 

for t E Z. 

Consider that the element at position in the matrix T(fc) is given by 

r(k) iJ =Cov(xt +k ,xi), 

thus modelling multivariate time series it holds that 

r(fc) = r(~kf, 

as matrices are not necessarily symmetric. This means that < Cov{X l t+k , X\ ) is 
in general not equal to <Cov{X 3 t+k , X'l) for k > 0. In the lead-lag setting the 
leading market's return 'yesterday' is strongly correlated to the lagging one's 
return 'today' but not vice-versa. 

Analogously to (1.13) the key to volatility contributions in this setting are 
the covariances of assets with the portfolio return. The following proposition 
gives the corresponding expressions for our setting. 

Proposition 2.6. Let (X t (X))tez denote the portfolio return when weighting 
the asset returns (Xt)t& by X, i.e. 

n 

x t {\) = Y,^xi forteZ, 

i=l 

then (X t (A))t e z is weakly stationary and it holds that 

7i(fc) := Cov(Xl +k ,X t (X)) = (r(fc)A) 4 , (2.8) 

for k = 0, 1, . . . and i = 1, . . . , n where (T(k)X)i denotes the the ith element of 
the vector T{k)X. 

Proof. The expectation of (X t (X)) t ^z is calculated straightforward. Further- 
more it holds that 

Cov(X t+k ,X t (X)) = Cov(X T X t+k ,X T X t ) = X T T(k)X, 

where T(k) denotes the covariance matrix of (Xt)tez for k = 0, 1, The above 

expression depends on the absolute value of the lag k only as 

A T r(-fc)A = A T r(fc) T A = A T r(/c)A, 
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which concludes the proof of weak stationarity. To prove (2.8) consider that by 
the bilinearity of covariance we get 



Cov(Xl +k ,X t (X)) = Cov(Xl +k ,X T X t ) 

n 

= Y,^ov{Xl +k ,X, 



for i = 1, . . . n and k = 0, 1, 



□ 



Using (2.8 1 we get useful expressions for the auto-covariance structure of the 
portfolio returns as well as risk contributions in this setting. 

Corollary 2.7 (Portfolio auto-covariance in the multivariate model). Let (X t (X))t^z 
be the portfolio return where the asset returns have covariance structure as given 
in (2.7) then the auto-covariance function (2.1) is given by 



7 (fc) = Cov(X T X t+k , X T X t ) = X T T{k)X = VK(*) 



(2.9) 



for i = 1, .. .,n and k = 0, 1, . . . where ji(k) is given in (2. 



Having calculated the auto-covariance function of (Xt(X))tez we find the 



volatility scaling for any d > 1 by Proposition 2.2 and Corollary [273] To conclude 
this section we analyze how to calculate contributions to volatility and derive 
formulas how these contributions change over time. 

Proposition 2.8 (Volatility contributions with serial correlations). Let (X t ) te % 
denote a weakly stationary process in W 1 and (X t (X))t^i the process of portfo- 
lio returns, then the volatility contributions by the Euler- allocation rule when 
holding this portfolio for d > 2 days are given by 



ai(X,d) 



X, 



d-l 



a(X,d) 



d li (0) + 2Y,(d-kH(k) , 



(2.10) 



k=l 



where 7i(fc) is given in ( 2.8 1 for i = 1, . . . , n 



fc = 0,l, 



The following corollary states how risk contributions scale in the above 
framework. 

Corollary 2.9 (Scaling of volatility contributions with serial correlations). Let 
<7j(A, d) as in Proposition 2.8, then for i = 1, . . . ,n 

(Ti{\,d) = <Ji(\)6(i,d), 



where a^(A) is the one-day volatility contribution as defined in (1.13) and the 
factor 5{i, d) is given by 



d-l 



8{i,d) 



7i(0) 



£(d-fc) 7i (fc) )/S(d), 



(2.11) 



k=l 



with ji(k) defined in (2.8) and 5(d) is the scaling factor for the respective port- 
folio volatility. 
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The proofs of the two statements above can be found in the appendix. 

Remark 2.10. Corollary |2.9| shows that in this modelling approach the relative 
volatility contribution changes depending on the holding period and the whole 
covariance structure of the returns involved. Thus assets whose relative risk 
contribution increases with the holding period can be identified. This is a feature 
that neither the model with uncorrelated asset returns nor the univariate model 



has. Finally, it is easily seen that (2.11) reduces to 



S(i, d) = —j= = Vd, for i = 1, . . . , n, 
v d 



in the case of no auto-correlations. 



Considering the results in Proposition |2.2| and 2J3 in general we have to 
estimate d — 1 auto-covariances respectively auto-covariance matrices to find a 
correct scaling for d days. Considering volatility per annum this corresponds to 
d — 1 = 249 or 251, which is clearly not appealing. In the following section we 
apply these findings to classical auto-regressive time series models which reduces 
the number of parameters tremendously. 



3 Scaling in ARMA and VARMA Models 

First we recall a few well-known definitions from the classical theory of time 



series analysis, see Brockwell and Davis 2002 , McNeil et al. 2005 , Taylor 



2008 . Note that for the ease of presentation we assume that all asset returns 
have zero mean. For daily returns this is usually assumed in risk management. 
We start with the basic building block in time series modelling. 

Definition 3.1 (White Noise). A process (Z t )t£Z is called (multivariate) white 
noise process WN(0, S) if it is covariance stationary and the covariance matrix 
function T(-) is given by 



T(k) = 



S, if k = 0, 
0, else, 



where £ is some positive-definite covariance matrix. 

By this definition there are neither serial cross-correlations between compo- 
nent series nor auto-correlations in the case of white noise. The only correlations 
exist at lag zero. Thus the uncorrelated portfolio case corresponds to a (multi- 
variate) white noise model. 

In order to analyze serial correlations we define Box- Jenkins models. 

Definition 3.2 (Box- Jenkins models). The weakly stationary process (X t )t^z 
with values in K™ is a called a VARMA (p, q) (vector auto-regressive moving 
average) process if it satisfies the following difference equation 

P Q 

k=l j=l 

where (Z t ) t6 z is WN(0, S) and fc = l,...,p and = 1, . . . , q are coef- 

ficient matrices in M. nxn . If Qj = for j = 1, . . . , q, then we denote such a 
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process by VAR(p), respectively VMA(g) if = for k = 1, . . . ,p . In the 
one-dimensional case such a process is denoted by ARMA(p, q) (auto-regressive 
moving average) and we write 



X t - 



p 

fe=i 



>kXt-k — 



9 



jZ t _j + Z t , for f g Z, 



(3.2) 



where {Zt)t&z is WN(0, <r), for some a > 0. 



In the following we will assume that the processes considered are causal. 
This means that they have a well defined representation as VMA(oo)-process, 
for a definition of causality we refer to the standard textbooks |Brockwell and| 



Davis p86 , Liitkepohl 2006 , Mills |1993|. 



Box- Jenkins models as defined above are frequently used to analyze financial 



time series and there exist various software packages (e.g. R Development Core 



Team 2010 , Matlab and mathematica) with functions to estimate the param- 
eters as well as to calculate the ACF ( |2.1[ ) respectively the ACF- matrices ( |2.7| . 

We quote the main propositions that link the parameters in (|3.1[) to the 
covariances given in ( |2.1[ ) respectively (2.7| (see for example [Brockwell and 
Davis[ |2002| Chapter 7] or |Liitkepohl| |2006 Chapter 2]). 



Proposition 3.3. The cross- covariances of a VARMA(p, q)-process, denoted by 



xz 



(k):=Cov(X u Z t _ k ), forkeZ, 



are given by 



and for k = 1, 2, 3 . 



T X z(k) = 0, fork<0, 
recursively by 



Txz(k)=J2^j T xz(k-j) + e k El k < q with T X z(0) = Z. (3.3) 

3=1 

After this preparation we can state the proposition that finally gives the link 



needed. See Mauricio 1995 for a proof and an algorithm for a fast and exact 



computation of the covariance matrices in the following proposition. 

Proposition 3.4 (Yule- Walker equations). The auto- covariance matrices T(k) 
for k — 0, 1, . . . of a causal VARMA(p, q) process are uniquely determined by the 
following system of linear equations 



r(k) - %r(fc - j) = 2 ©3-rxz(j - k). 

3=1 3=fc 



(3.4) 
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In the following we derive concrete formulas for scaling in the univariate 



case. 



Example 2 (MA(g)-model). For the MA(g)-model the formula for ACF is ex- 
plicitly found in many textbooks (e.g. |McNeil et al. 2005| Chapter 4.2]) fur- 
thermore Proposition |3 . 4| gives , setting 9q = 1, 



7 (0) =<r 2 X>J, and, 

3=0 
9-1*1 

1 {k)=a 2 J2 e A+\k\, for |*| = 1, 2,..., g. 

3=0 



(3.5) 



Obviously j(k) = for |fc| > q. 



We can plug the acf (3.5| into the formula (2.5) and see that the correct 



scaling factor in this setting is given by 



\ 



Min(d-l,q) 

d + 2 (d-k) 

k=l 



'jVj+\k\ 



ELo 6 3 



Example 3 (AR(l)-model). Incorporating an auto-regressive coefficient we con- 
sider an AR(l)-model whose ACF^Jis given by McNeil et al. 2005 Chapter 4.2] 



7W = P 



1 - ■ 



for Ifel = 0,1,2,... 



Note that in this case 7(fc) ^ for k 6 Z and we have to apply standard 
summation calculus to get the following scaling factor 



(5(d) = Jd + 2- 



r (d(l + 



3.1 Closing time problem 

Before we consider the closing time problem it is essential to understand the in- 
terplay between multivariate time series and portfolios constructed from these. 
Creating a portfolio return (A t ) te z from asset returns modelled by a multivari- 
ate time series (X t ) teZ using the weights A by X t — X T X t means to apply a linear 
transformation to X t . To understand this transformation we quote the following 



proposition on VMA(g)— processes Lutkepohl 2006 Proposition 11.1] 



Proposition 3.5 (Linear transformation of a VMA(g)-process). Let (Z t ) te z be 
an n-dimensional white noise process with nonsingular covariance matrix £ and 
let 

q 

X t =J2 e 3 Z t-j+Z u forteZ, 



3 With increasing order of the models handy formulas for the auto-covariance function 
can not be provided any more, but the function ARMAacf in the software package R (see |R| 



[Development Core Team] |2010| ) or similar implementations may be used to find the ACF 
after having estimated the ARMA-coefncients. 
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be an n- dimensional invertible VMA(q) -process. Furthermore, let F be an (M x 
n) matrix of rank M. Then the M- dimensional process 3^ = FX t has an invert- 
ible VMA(q) representation, 



y, 



q 



BjZt-j + Zt, fort el, 



where (Z t )tez is M-dimensional white noise with nonsingular covariance matrix, 
the Qj are coefficient matrices and q < q. 

This proposition can be applied to analyze the situation if we form a portfo- 
lio. Then F = X T is 1 x n and we get a moving average process of order equal or 
less the order of the vector process. For the closing time problem q = 1 which 
leads to q = 1 and we can summarize our findings for this problem so far: 

Observation 1. The multivariate process of closing-time returns of assets 
traded in different time zones can be modelled as a VMA(l) process of the form 

x t = e 1 z t _ 1 + z u forteZ. 

Creating a portfolio with asset weights X results in a MA(1) process of the from 
X T X t =: X t = 0i.Zi-i + Z t , for t <= Z. 
For a MA(l)-process and the scaling constant of Example [2] simplifies to 



5{d) = Jd + 2(d- 1) 



1 



for d > 1. 



(3.6) 



This allows a top-down modelling of the returns if we think that the only 
auto-correlations in the portfolio come from the closing time problem. We 



just have to fit the parameters 9 and a by (3.5) and apply (3.6) for correct 



scaling of portfolio volatility. Alternatively, we can of course directly estimate 



the autocorrelation p(l) of portfolio returns and plug it into (2.5). However, 
estimating the full model would give as more insight. 



The VMA(l)-model In this paragraph we derive the details of the formulas 
of Section[2]for the VMA(l) -case. Below in Example [4] we apply these formulas 
to the closing time problem started in Example [l] By ( |3.3[ ) and (3.4) we get 
the following for the VMA(l)-model: 



which gives 



r(o) ^e^e^ + s, 
r(i) = e 1 s, 

7 (0) = A T (GiSef + E) A, 
7 (1) = A T (6 1 S)A ! 



and 



and 



(3.7) 



(3.8) 



for the portfolio time series by (2.9). Thus the scaling constant (2.5) is given by 



5{d) = Jd + 2(d- 1) 



A T 9i£A 



A T (0iE6f + E) A' 



(3.9) 
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and the scaling of contributions (2.11 ) is given by 



5{i,d) = d + 2(d- 1) 



(gigA)j 

((e^ef + e)a). 



/5(d), 



(3.10) 



where 6(d) is calculated in (3.9) above. Again as in Example [2] one can estimate 
r(0) and r(l) and plug them into (2.11 1. But the estimator for Oi and £ 
or transformations of it will reveal interesting structures (see e.g. Tsay 2005 



Chapter 8.2.1 Reduced and Structural Forms] or Liitkepohl 2006 Chapter 2.3.2 
Impulse Response Analysis]). 



Scaling volatility of closing-time returns in VMA(l) compared to scal- 
ing volatility of contemporaneous returns In this paragraph we have a 
closer look at the Newey-West estimator (1.9) and the naive estimator (1.8) in 



the the context of volatility scaling. Having estimated the lag-zero covariance 
matrix of asset-returns T(0) and the lag-one covariance matrix T(l) we find the 
volatility of the d days portfolio closing-time return by (2.9) as 



<r(\, d) = ^dX T T(0)X + 2(d - l)A T r(l)A, for d > 1. (3.11) 

Considering the contemporaneous returns together with their covariance ma- 

and assuming zero autocorrela- 



trix given by the naive estimator £ from ( 1 



tions among contemporaneous returns the d days volatility of the contempora- 
neous portfolio return cx(A, d) is given by 

cr(A,d) = VdX T tiX 



^/dA T (r(o) + r(i) + r(i) T )A 



= JdX T T(0)X + 2dA T r(l)A, for d > 1. 



(3.12) 



Considering the ratio of the scaled volatility of the portfolio closing-time re- 



turn (3.11) over the scaled volatility of the contemporaneous portfolio return (3.12) 



we see that this quantity converges to 1: 



Jim w^'ll = Jim 

d— s-oo <7^A, d) d— >co 



y/ dA T r(0)A + 2(d - l)A T r(I)A 
v/dA T r(0)A + 2dA T r(l)A 



= 1. 



(3.13) 



Thus for large d the risk figures for the two procedures coincide. 



However, applying the Newey-West estimator up to lag 1 (1.9) the limit of 
the corresponding ratio is given by 



lim 

d— ¥oo 



V / dA T r(0)A + 2(d - l)A T r(T)A 

yj dA T r(0)A + dA T r(l)A 
VA T r(0)A + 2A T r(l)A 
v /A 7 T(0)A + A T r(l)A 



^1. 



Thus the volatility estimates by the Newey-West estimator (1.9) for contem 



poraneous returns, in this set-up, do not coincide with the result of applying 
the VMA(l)-model for closing-time returns nor with the result of applying the 
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naive (but in this case correct and useful) estimator (1.8). This shows that the 
compatibility of long-term risk estimates in the two notions, closing-time return 
and contemporaneous return, depends on the estimator used for the covariancc 
matrix of contemporaneous returns. 



As already mentioned the estimator £ defined in (1.8) can, in general, be 



invalid which by (3.13) questions the VMA(l)-model. In this case we should 
analyze the situation in more detail and make sure that the assumption that the 
closing-time problem is the only source of auto-correlation is acceptable. In the 
following we apply our findings to the data from Example [I] and all estimators 
considered are mathematically valid. 

Example 4 (Example [I] continued). The conclusion of Example [l] is that we can 
model the closing-time returns as VMA(l)-process. We apply this to the data 



of Example [TJ Note that we can not solve (3.7) for Oi or £. Thus we apply 



maximum likelihood estimation for this task as it is provided in the R-package 
DSE |Gilbert| [2006 or later] . Note that positive definiteness of the covariance 
matrix of residuals of the VMA(l)-model is assured in the estimation procedure 
applied This fact is required for valid estimators in (3.7). 



On the other hand we can directly estimate T(0) and T(l) from the data, 
knowing that T(k) — for k > 2. But looking at 0i and £ gives us some 
insight in the problem at hand. As risk and risk contributions is the focus of 
this article we do not go through the whole impulse-response analysis but refer 
to Tsay 2005 Chapter 8.2.1 Reduced and Structural Forms] or Liitkepohl 



2006 Chapter 2.3.2 Impulse Response Analysis]) or the generalized impulse- 



response of Pesaran and Shin 1998] . 

In Figurc[3]we see the scaling constant for the portfolio time-series by the just 
mentioned MLE estimation of the VMA(l)-model of the seven assets, a MA(1)- 
model estimated directly on the portfolio time-series and the square-root-rule. 
The reason why the upper lines do not match perfectly are estimation errors 
but this plot gives us some confidence for the MLE estimators for Oi and £. 

Considering Table[T]we see the contributions to volatility p. a. on one hand by 
assuming zero serial correlations between the assets in the portfolio and on the 
other hand by modelling the asset returns in the portfolio as VMA(l)-process. 
Applying the square-root-of-time rule to this portfolio we get a volatility p. a. of 
20.07% while the volatility p. a. increases by approximately 17% to 23.46% in 
the VMA(l)-model. Concerning the analysis of sources of risk note that the risk 
contribution by the geographically most distant market Japan of 1.58% looks 
quite small when ignoring serial cross-correlations but it increases to 3.84% 
taking them into account. The risk contributions of the markets leading the 
portfolio do not change dramatically. The increase of the portfolio volatility can 
be attributed to the Asian assets whose risk contribution increases significantly 
as serial cross correlations to the US and Europe are taken into account. This 
example shows that not only the accuracy of total volatility of the portfolio 
increases but also the attribution to single assets fits economic considerations 
much better! 



We acknowledge personal communications with the author, Paul Gilbert, on this topic. 
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50 100 150 200 250 

d 



Figure 3: Scaling constants 5(d) for ranging d for a full VMA(l)-model (solid, 
gray), a univariate MA(l)-model (dashed, black) and the square-root-of-time 
(dotted, black). 



Asset 


Currency 


Exposure 


Square-root rule 


VMA(l) 


Difference 


Portfolio 


EUR 


105% 


20.07% 


23.46% 


3.38% 


Topix 


JPY 


15% 


1.58% 


3.84% 


2.26% 


H-shares 


HKD 


15% 


3.56% 


5.40% 


1.84% 


DJ Euro Stoxx 50 


EUR 


15% 


3.28% 


2.96% 


-0.32% 


Swiss Market 


CHF 


15% 


2.24% 


2.11% 


-0.13% 


JSE TOP 40 


ZAR 


15% 


2.62% 


2.96% 


0.34% 


Russell 2000 


USD 


15% 


3.94% 


3.43% 


-0.51% 


NASDAQ 100 


USD 


15% 


2.85% 


2.75% 


-0.11% 



Table 1: Contributions to volatility p. a. in a the global portfolio of Example [T] 
applying the square-root-of-time rule and modelling a VMA(l)-process. 



3.2 Genuine auto-correlations 

We conclude our theoretical study of auto-correlated portfolio returns by a 
short detour to genuine auto-correlations. In the literature studies can be found 
(see e.g Anderson et al. 2005 and references therein) which provide evidence 



that, besides the spurious effects of non-contemporaneous trading, genuine ef- 
fects such as partial price adjustment and time- varying risk premia can lead 
to genuine auto-correlations in asset returns. We stay in the class of vector- 
autoregressive models but focus on VAR(p), especially, VAR(l)-models in this 
context as opposed to the VMA(l)-model for the closing time problem above. 

Before we consider an example of first order genuine auto-correlations we 
go one step deeper into understanding the interplay between multivariate time 
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series and portfolios constructed from these. An important result on linear 



transformations of VARMA(p, q) processes is the following Lutkepohl 2006 
Corollary 6.1.1]: 

Theorem 3.6 (Linear transformations of VARMA(p, q) processes). Let (X t )tez 
be an n- dimensional, stable, invertible VARMA(p,q) process and let F be an 
M x n matrix of rank M. Then the process (FX t ) te % has a VARMA(p,q) 
representation with 

p < np and q < (n — l)p + q. 

The above theorem tells us that a portfolio which is a simple linear transfor- 
mation of a VARMA(p, q) process can not be guaranteed to have an ARMA(p, q) 
representation of the same order. Furthermore it is important to note that 
as a consequence the class of VAR(p)-models is not closed with respect to 
linear transformation as, in general, the result of the transformation can be 
some VARMA(p, q) process with q > 0. Concerning multivariate models we 
nevertheless focus on VAR(p)-models due to known identification problems of 
VARMA(p, q)-models if q > (see for example |Lutkepohl| [2004] ) . 



As a preparation for our key result on portfolios constructed out of VAR(p)- 
process we state the following proposition: 

Proposition 3.7 (AR portfolios built from VAR processes). Let (X t )tez be a 
VAR(p) -process in 1™ for n > of the form 

p 

and let A = (Ax, . . . , A„) be a vector of weights in K™. Then the portfolio process 
(A t (A)) te z is an AR{p)-process of the form 

p 

X t (\)^J2^X t -k(X) + X T Z t (3.14) 



fc=i 



if and only if 



A T $fe = fc A T fork = l,...,p. (3.15) 



Note that as in Corollary 2.7 (A Z t ) t€ % is clearly a white noise process. The 



full proof of the proposition can be found in the appendix. Condition (3.15) 



means that only portfolios with A being an eigenvector of all coefficient matrices 



of the VAR(p)-process admit the AR(p) representation (3.14| which one could 
expect to hold in general, at first glance. The following corollary concludes these 
considerations. 

Corollary 3.8 (AR portfolios built from VAR processes). In the setting of 



Proposition 3.7 the process (X t (X))t<£Z is an AR(p)-process for any portfolio 
weighting A € R" if and only if the coefficient matrices of the VAR-process are 
diagonal and of the following form 

$fc = (j>kl n , fork=l,...,p. 
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The consequence of the above corollary is that modelling genuinely auto- 
correlated assets we will in general not observe portfolio returns consistent with 
an AR(l)-model. This would only be possible if the coefficient matrix $1 were 
of the form 

(a ... 0\ 
a ... 

: '•• : 
\0 a) 

for some fixed value for a - all the same for each asset. The conclusion is that 
the weighted VAR(l)-model is richer than an AR(l)-model. 



The VAR(l)-model After these general considerations we focus on the VAR(l)- 
modcl of the form 

Xt = $iXt-i + Zt, forieZ, 

since this model will be the natural choice to capture genuine serial correlations. 
In this case we get the following expressions for the covariance matrices 



which gives 



r(0) = S(J - Qjy 1 and 

r(jfe) = $*r(o), for k > i, 

7 (0) = A T r(0)A, and 

j(k) = A T (#Jr(0)) A, for k > 1, 



(3.16) 



(3.17) 



for the portfolio time series by ( |2.9[ ). Using (3.161 and (3.17) the scaling con- 
stant (2.5| for d > 1 is given by 

S(d) 



\ 



d-l 



k=l 



A T $jT(0)A 
A T T(0)A ' 



(3.18) 



and the scaling of contributions (2.11 1 is given by 

d-l 



6(i,d) 



1 



5(d) 



d + 2^2(d-k) 



fc=i 



(T(0)A)i 



(3.19) 



where 5(d) is calculated in (3.18) above. In contrast to the closing-time problem, 
in this case, we would have to estimate d—1 lagged covariance matrices T(fc), k = 
1, ... ,d — 1 if we wanted to plug them into (2.11 ) directly. This is clearly not 
feasible for large values of d which justifies the use of a specific time-series model 
in these cases. 

The following concrete example illustrates the above issues. We further- 
more analyse the trade-off when approximating such a portfolio with first order 
genuine auto-correlations with an AR(l)-model, although it is not theoretically 
justified. 

Example 5. Consider a portfolio consisting of two contemporaneously traded 
assets A and B with annual volatilities of 25% and 20% and a correlation of 
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70%. Asset A has a negative genuine first order autocorrelation of —5% while 
asset B exhibits 2.5% genuine first order autocorrelation. Note that these values 



are consistent with findings in Anderson et al. 2005 . Thus we consider the 
following covariance matrices 

T(0) = diag (0.25 0.2) Q ? °A diag (0.25 0.2) ^ and 



r(l) = diag (0.25 0.2) 



0.05 \ ,. / no , no x 1 
diag (0.25 0.2)—. 







0.025 



We can model these two assets by a VAR(l)-model and calculate the coefficient 
matrix $i by (3.16) and get 



$i = r(i)r(o) _1 = 

Modelling a portfolio with the weighting A 



0.0980 0.0858\ 
0.0275 0.0490 J 



(§, \) T we expect an ARMA(2,1) 
process by Theorem 3.6 However, for a pure AR(l)-model the coefficient 4>i is 
given by 

&= 7 (l)/ 7 (0) = ^!= -0.0123. 

Using the results for scaling in univariate models from Example [3] and Equa- 
tion (3.18) we compare the resulting scaling constants in Table [2] We see that 
the AR(l)-model performs well especially for shorter holding periods in approx- 
imating the result of the multivariate model. However, such a simple approxi- 
mation should be used with care. Furthermore the VAR(l)-model tells us more 
details about the risk contributions as wc will sec below. 



d 


VAR(l) 


AR(1) 


SRTR 


2 


1.405 


1.405 


1.414 


5 


2.218 


2.214 


2.236 


10 


3.134 


3.127 


3.162 


30 


5.427 


5.412 


5.477 


90 


9.398 


9.372 


9.487 


20 


15.662 


15.619 


15.811 



Table 2: Volatility scaling factors 8(d) for VAR(l), AR(1) and the square-root 
of time for various holding periods d. 



We conclude this detour on genuine auto-correlations by an analysis of the 
relative risk contributions, i.e. risk contributions in percentage of total volatility. 
Applying (3.19) to our example we see in Table[3]that the contribution of asset 
A is dominant as it has the higher volatility. But with increasing holding period 
the relative risk contribution of asset A decreases which reflects its negative 
auto-correlation and the positive auto-correlation of asset B. This is a feature 
that only the multivariate approach can offer. 



4 Conclusions 



In this article we first clarify the notion of closing-time returns and contempo- 
raneous returns in global portfolios. In Example [I] we illustrate these notions 
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d 


a(d,A) 


a(d.B) 




a{d) 


a(d) 


1 


56.52 


43.48 


2 


55.39 


44.61 


5 


54.77 


45.23 


10 


54.56 


45.44 


30 


54.42 


45.58 


90 


54.37 


45.63 


250 


54.36 


45.64 



Table 3: Relative risk contributions (percentage) of assets A and B for various 
holding periods d in the VAR(l)-model. 



in a setting of time-shifted multivariate Brownian motion. Serial correlations 
naturally occur when analyzing portfolios of geographically diversified assets 
traded in distant time zones and we motivate the use of a VMA(l)-model for 
the closing-time returns. We then address the problem of calculating portfolio 
volatility of closing-time returns for holding periods of more than one day. 

We show that ignoring serial correlations leads on one hand to biased esti- 
mates of volatility and on the other hand to misleading risk contributions as 
Example [4] illustrates. 

We propose formulas for calculating accurate volatility scaling modelling the 
portfolio closing-time return as a univariate process as well as in a multivariate 
setting. Moreover in the multivariate setting we also provide explicit formulas 
for genuine risk contributions that take the time series structure of the assets 
involved into account. Modelling the asset returns as a vector moving average 
process of order one we derive handy formulas and perform a complete analysis 
of risk and risk contributions and compare this approach to the Newey-West 
estimator of contemporaneous returns and another simple but useful estimator 
in the same spirit. 

Finally we take a short detour to genuine auto-correlations and propose the 
application of a VAR(l)-model to tackle this problem. 

Applying the findings of this article to the calculation of the tracking error, 
i.e. the volatility of the additional return of the portfolio above a given bench- 
mark, can improve the analysis of relative risk which is often an aim in asset 
management. 

Besides the analysis of market risk our findings can be applied to portfolio 
optimization as well as portfolio construction techniques such as risk-parity (also 
known as equally- weighted risk contributions, see Maillard et al. |2008| ) where 



risk contributions by assets are the driving input. As another direction of further 
research the findings of this article may also be applied to the VEC specification 
of multivariate GARCH models, since they admit a VARMA representation 
(see 



Lutkepohl 2006 ) 



Appendix: Proofs 



Proof of Proposition\2. 2\ and Corollary 2.3 Considering that er(A, d) is the square- 



root of Var(y^ _ 1 Aj(A)) and writing down this variance in matrix form using 
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weak stationarity we get 



i=l 



/Cov(X 1 {X),X 1 {X)) 
Cov(X 2 (X),X 1 {X)) 



VCou(X d (A),X 1 (A)) 

/ 7(0) 7(1) 
7 (1) 7 (0) 

\j(d-l) 7(d"2) 



Co«(X 1 (A),X d (A))\ 
Cou(X 2 (A),X d (A)) 

C™(X d (A),Jf d (A))7 

7 (d-l)\ 
7 (d-2) 

7(o) y 



where 7(-) denotes the auto-covariance function of {X t (X)) t £z and 1 = (1, . . . , 1) T 
Summing up along the diagonals and using symmetry we get Equation (2.3 1. 
For proving (2.5| note that the auto-covariances can be expressed in terms of 

n; 

^ 1 (k)= 1 (0)p(k) = cr(X) 2 p(k), 



the auto-correlation and the variance in the following sense 
7(fc) 



p(k) 



7(0) 



for k = 0, 1, . 



□ 



Proof of Proposition \2.£\ and Corollary \2.£\ Following the Euler allocation rule 
the volatility contributions are given by 



<Ti(\,d) = Aj 



ct(A, d) 



where we can calculate the denominator <r(A, d) by (2.3 1 and (2.9). For the 
numerator we get 



CovC£xiY,x k (\)) 



fc=l 



k=l 



/CoviXiX^X)) 
Cov(XiXi(X)) 



yCoviX^X^X)) ... 

/ 7,(0) 7,(1) 
7,(1) 7,(0) 

\-Yi(d-l) n{d-2) 



Cov{X{,X d (X))\ 
Cov{XiX d (X)) 

Cov(X d ,X d (X))J 

.. 7 ,(rf-l)\ 
.. Hid -2) 

■■ 7,(0) J 



where 1 = (!,...,!) and the 7 j(fc) are given in (2.8) for i = 1, . 



.,n 



and 



fc = 0, 1, . . . Now again use the symmetries and sum up along the diagonals to get 
the result. To prove the corollary, recall that a(X, d) = a(X)S(d). Plugging this 
into the denominator of (2.10) and recalling that the one day risk contribution 
is given by <7j(A) = A, we get 



cr,(A, d) = A. 



7i(0) 

(7(A) 



7,(0) 



5>-A) 7< (fc) /*(d), 



fe=i 



which gives the form of the factor S(i,d) for i = 1, 



□ 
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Proof of Proposition\3.7\ and Corollary \3.S\ Considering Proposition |3.7| we use 



that by assumption (A t (A))t e z is an AR(p)-process as in (3.14) and the identity 
A At = Xt'. 

\ T X t = X t (\) 
p p 

>?<$>kXt-k + X T Z t = J2 <i>k\ T Xt-k + \ T Z t 

k=l k=l 
P 

«■ A T ^($ fc -0 fc /„)A' t _ fe = 0. 
fe=i 

As (X t )t£Z takes values in R n this forces 

A T ($ fe -<Mn)=0 forfc = l,...,p. 



To prove Corollary |3.8| consider that the above condition is true for arbitrary 
A G W 1 if and only if the rank of the matrix — is zero for k = 1, . . . ,p 
which concludes the proof. □ 
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